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Abstract. The time evolution of pseudorapidity distributions of produced charged hadrons in d+Au col- 
lisions at -y/sjvjv = 200 GeV is investigated. Results of a nonequilibrium-statistical Relativistic Diffusion 
Model with three sources are compared with a macroscopic "bounce back" model that does not allow for 
statistical equilibration at large times, but instead leads to motion reversal. When compared to the data, 
the results of the diffusion approach are more precise, thus emphasizing that the system is observed to be 
on its way to thermal equilibrium. 

PACS. 25.75.-q Relativistic heavy-ion collisions - 24.60.Ky Fluctuation phenomena 



1 Introduction 

It has recently been shown that rapidity distributions of 
net protons [TJ[2] and pseudorapidity distributions of pro- 
duced charged hadrons [3 ,4, 7] in relativistic heavy-ion col- 
lisions can be described very precisely in a nonequilibrium- 
statistical Relativistic Diffusion Model (RDM) with three 
sources. There are two sources with initial rapidities close 
to the beam values that arise mainly from the valence 
quarks, and a third midrapidity source that is mostly due 
to gluon-gluon collisions, which reaches statistical equi- 
librium with respect to the variable pseudorapidity dur- 
ing the time evolution in the diffusion model. The parti- 
cles in this equilibrium source move collectively with very 
large velocities [2] close to light velocity similar to a blast 
wave They may have hadronized from an equilibrated 
quark-gluon plasma. 

In case of the asymmetric system d+Au at the highest 
RHIC energy ^snn — 200 GeV, we could demonstrate 
that this analytical model of a many-particle system on 
its way to statistical equilibrium reproduces the compli- 
cated features of the experimental distribution functions 
to a high degree of accuray. Signatures such as the steeper 
slope of the distributions in the deuteron direction as com- 
pared to the gold direction, the gradual change of the 
pseudorapidity distributions with ccntrality, and the dom- 
inance of particle production in the Au-like rapidity region 
towards more central collisions arises in a natural way 
from the nonequilibrium-statistical description. The mea- 
sured pseudorapidity distributions remain rather asym- 
metric because strong interaction stops before the system 
reaches statistical equilibrium. 



In this work we compare the time evolution of the pseu- 
dorapidity distribution in the nonequilibrium-statistical 
RDM [1,2,3,1] with a schematic macroscopic model that 
does not incorporate statistical equilibrium as a limit for 
large times (except for the midrapidity source), but in- 
stead leads to a separation of the beam-like distribution 
functions in pseudorapidity space. In the course of their 
time evolution, the mean values of the beam-like distri- 
bution functions change their signs. This corresponds to 
motion reversal, or a " bounce back" of the produced par- 
ticles in the beam-like sources with respect to the beam 
directions. Such a model does not lead to statistical equi- 
librium in the system for large times since the partial dis- 
tributions separate in pseudorapidity space. 

In case of an asymmetric system such as d+Au investi- 
gated here, it is possible to decide whether agreement with 
the data occurs before or after motion reversal. Although 
a corresponding distinction is not possible for symmetric 
systems such as Au+Au, the results for d+Au will indicate 
that symmetric systems behave in an analogous way. 

The Relativistic Diffusion Model is reconsidered briefly 
in the next section. Afterwards, the bounce back approach 
is outlined, and both models with their corresponding time 
evolution are then compared to the available PHOBOS 
minimum-bias [9 and centrality dependent [TD] d+Au data 
at ^/s~mn — 200 GeV, and to each other. Such a compar- 
ison is conceptually quite important because a preference 
for the diffusion case would emphasize that the relativistic 
systems are indeed observed to be on their way to thermal 
equilibrium. 
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2 Nonequilibrium-Statistical Approach 

The nonequilibrium-statistical approach is based on diffu- 
sion equations which eventually lead to statistical equilib- 
rium for large times. In the three-sources model, we have 
successfully used linear Fokker-Planck equations (FPEs) 
HJEE] for the components Rk(y,t) of the distribution 
function for produced charged hadrons in rapidity space 



dt 



R k (y,t) 



id 

T y dy 



(y-yeq)-Rk(y,t) 



dy 2 



D k v -Rk(y,t) 



(1) 



with the rapidity y = 0.5-ln((£ , +p) /(E—p)). The diagonal 
components D y of the rapidity diffusion tensor contain the 
microscopic physics in the respective beam-like (k = 1,2) 
and central (k = 3) regions. They account for the broaden- 
ing of the distribution functions through interactions and 
particle creations. The off-diagonal terms of the diffusion 
tensor account for correlations between the three sources. 
They are expected to be small, and we neglect them in 
this work. The rapidity relaxation time t v determines the 
speed of the statistical equilibration in y-space. 

Whereas an equilibrium-statistical description of rel- 
ativistic heavy-ion collisions works well for certain ob- 
servables such as abundance ratios of produced hadrons 
(thermal model) , a nonequilibrium-statistical formulation 
as outlined here is required for distribution functions and 
other more sophisticated observables. A description of the 
dynamics should be part of such a formulation. In the 
present work, the dynamics enters through the time vari- 
able in the FPE that describes the evolution of the three 
sources for particle production in rapidity space. 

For time to infinity, both the mean values and the vari- 
ances of the three sources reach the equilibrium values 
with respect to the variable rapidity. Then the incoherent 
sum of the three distributions represents the overall statis- 
tical equilibrium distribution in rapidity space and after 
Jacobi-transformation, in pseudorapidity space. Compar- 
ison with the data for d+Au at the highest RHIC energy 
shows, however, that strong interaction stops before this 
equilibrium distribution is reached. 

In the nonequilibrium-statistical model as well as in 
the bounce back approach, the initial particle production 
in the three sources is assumed to occur with full strength 
at very short times. In the diffusion case, we use in this 
work the initial conditions Ri^iy^t = 0) = 6(y ± y ma x) 
with the maximum rapidity y max = 5.36 at the highest 
RHIC energy of ^Jsnn — 200 GeV (beam rapidities are 
2/1,2 = TVinax), and R 3 (y,t = 0) = S(y). A midrapid- 
ity gluon-dominated symmetric source had also been pro- 
posed by Bialas and Czyz [6j. It may originate from a 
thermally equilibrated quark-gluon plasma. 

These initial conditions are slightly modified for the 
central source as compared to our previous investigation 
[4], where we used i?3(y, t = 0) = 6(y — y eq )- The motiva- 
tion is twofold: We want to use identical initial conditions 
when comparing diffusion and bounce back approach, and 
we want to investigate whether the modified initial condi- 



tion allows to maintain the quality of the agreement with 
the data in the diffusion approach. 

The new initial condition for the midrapidity source 
corresponds to initial (t=0) particle production at rest, in- 
dependently of the mass of the collision partners: the third 
source is created at y=0. Hence, this gluon-dominated 
source is initially insensitive to the mass distribution in 
the system, and it is only at larger times - when the par- 
ticles in this source have already been created - that the 
drift towards the equilibrium value sets in. Since there 
are very few valence quarks in the midrapidity region at 
the highest RHIC energy, this initial condition is probably 
more realistic than the one we used in [3]. 

The mean values in the three sources have the time 
dependence 

<ydiA t ) >= yeq[l-exp(-t/Ty)}^y max exp(~t/T y ) (2) 



for the sources (1) and (2), and 



< yd3(t) >= 2/ eg [l - exp(-t/r y )} 



(3) 



for the moving central source. The three mean values reach 
the equilibrium value for time to infinity. In our previous 
RDM-calculation [4 with slightly different initial condi- 
tion, the mean value of the central source was at the equi- 
librium value < j/3 (£) >= y eq independently of time, thus 
assuming instant equilibration in this source regarding the 
mean values. The variances are obtained as in [5]. It will 
turn out that for d+Au at the highest RHIC energy, we 
can not determine from a comparison with the data which 
of the two possibilities for the initial conditions of the 
central source is more realistic because the \ 2 is nearly 
identical in both cases. 

We follow the subsequent diffusion-model time evolu- 
tion in pseudorapidity space up to the interaction time 
T int , when the produced charged hadrons cease to inter- 
act strongly. The quotient Ti n t/T y is determined from the 
minimum x 2 with respect to the data, simultaneously with 
the minimization of the other free parameters - in partic- 
ular, the variances of the three partial distribution func- 
tions, and the number of particles produced in the cen- 
tral source. In the nonequilibrium-statistical approach, the 
equilibrium value of the rapidity and its dependence on 
centrality is calculated from energy and momentum con- 
servation in the system of participants as 

Ve q (b) = 

1 ln < mj(b) > exp(-y max )+ < m%(b) > exp(y max ) 

2 < m%(b) > exp(-y max )+ < mf(b) > exp(y max ) 

with the transverse masses < mf 2 (b) >= \f[m\ 2 (b) + 

< pr > 2 ), and masses mi j2 (6) of the "target" (Au)- and 
" projectile" (d)-participants that depend on the impact 
parameter b. The average numbers of participants 

< Ni^(b) > from the Glauber calculations reported in [ID] 
for minimum bias d + Au at the highest RHIC energy are 

< Ni >=6.6, < N 2 >=1.7, which we had also used in [U 
|S]. With < p T >= 0.4 GeV/c the result is y eq = -0.664, 



with < pt >= 1 GeV/c we obtain y eq 



-0.60. 
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The average numbers of charged particles in the target- 
and projectile-like regions are proportional to the 

respective numbers of participants N\^, 



Ni. 



{N tot 



(JVi + N 2 ) 



(5) 



with the constraint 



N 2 ch 

riot 



N2- Here the 



total number of charged particles iV'£ is determined from 
the data. The average number of charged particles in the 
equilibrium source is a free parameter that is opti- 
mized together with the variances and Ti n t/T y in a x 2- fit 
of the data using the CERN minuit-code. Due to the ac- 
curacy of the data, these five free parameters (for both 
scenarios, diffusion and bounce back) are determined with 
great precision. 

The FPE is solved analytically as outlined in [U[S] , the 
result is converted to pseudorapidity space, and compared 
to data, figjl] right-hand column. Here the time evolution 
parameter p in the numerical calculation is defined as Q 

p= 1 - exp(-i/T„). (6) 

We have shown in [H|5] that with this RDM approach, 
the centrality dependence of the measured pseudorapidity 
distributions [10] from central to very peripheral collisions 
can be modeled in considerable detail. For peripheral col- 
lisions, the asymmetry of the overall distribution is not 
yet pronounced because here the d- and the Au-like par- 
tial distributions are similar in size due to the small num- 
ber of participants. Towards more central collisions, the 
number of gold participants rises, and the corresponding 
partial distribution of produced particles becomes more 
important. In addition, the distributions drift towards the 
equilibrium value. Both effects produce the asymmetric 
shape, which is also seen in minimum-bias. 

The minimum-bias result that we present here in more 
detail also shows the asymmetric shape, which is very well 
reproduced in the diffusion calculation. At larger values 
of the time evolution parameter p, all three subdistribu- 
tions tend to become symmetric in y with respect to the 
equilibrium value y eq , indicating the approach to thermal 
equilibrium. At p=0.999, the equilibrium state is already 
closely approached. The slight asymmetry is due to the 
conversion from rapidity- to pseudorapidity space which 
tends to produce a dip at i] = 0. For time to infinity, 
statistical equilibrium in pseudorapidity space would be 
reached. 

It is interesting to determine actual values for the diag- 
onal components of the diffusion tensor from the expres- 
sion for the variances, eq.(lll, provided the interaction 



time is known. As an example, for an interaction time of 
7fm/c in d+Au and T int /r a =0.78, the resulting values of 
the three diffusion coefficients in minimum-bias are Dy' 2 ' 3 
= 1.55, 0.46 and 0.98 c/fm. These are effective values as 
explained in ref.[2] because they include the influence of 
collective expansion in addition to the statistical fluctua- 
tions. 



There is a difference of a factor of two in the exponent 
as compared to the definition of p used in [S], which causes 
different t/r y values for given p. 



Table 1. Produced charged hadrons in minimum-bias d + Au 
collisions at y'sjvjv = 200 GeV, yi.2 = T 5.36 in the "bounce 
back" model (1st line), and in the Relativistic Diffusion Model 
(2nd line). The equilibrium value of the rapidity in the RDM 
is j/ e( j, the time parameter (see text) is p, the corresponding 
value of interaction time over relaxation time is Ti n t/T v , the 
variance of the central source in y— space is erf. The number of 
produced charged particles is for the sources 1 and 2 and 
N^ h for the central source, the percentage of charged particles 
produced in the midrapidity source is n^. h . 



3 Bounce Back Approach 

In this schematic model of motion reversal in pseudora- 
pidity space, the beam-like partial distribution functions 
change their signs in the course of the time evolution of 
the collision. The behaviour of the system is approximated 
through a linear partial differential equation which is sim- 
ilar to a Fokker-Planck equation, but differs in that it 
does not lead to statistical equilibrium in the system for 
large times. Instead it causes a separation of the beam-like 
sources in rapidity space: 
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—R k (y,t) = — — 
at t v dy 



(y-yk)-Rk(y,t) 



dy 2 



D k v -R k {y,t) 



(7) 



Here the drift term does not account for the equilibration 
towards y e q(b) as in the FPE, but rather describes the 
motion of both peripheral distributions in rapidity space 
towards smaller absolute values of the rapidity, and even- 
tually into the other hemisphere of rapidity space - which 
means motion reversal, or bounce back. The drift term 
that induces this behaviour of the beam-like partial dis- 
tributions is {y - y k ) = (y T Umax) for k=l,2. The drift 
term for the midrapidity source is (y — y^) = y. 

We use the same t v as in the RDM case in order to have 
a direct comparison of the two scenarios. In the bounce- 
back scenario, there is no relaxation towards statistical 
equilibrium and hence, r y should not be considered to be 
a relaxation time in this case, it is just a time parameter 
that controls the speed of the movement in y- (or ?y-)space. 

Because the three-component system does not approach 
statistical equilibrium with respect to the variable rapid- 
ity for time to infinity, the underlying system of partial 
differential equations for the bounce back case differs in 
an important aspect from the ordinary Fokker-Planck, or 
Uhlenbeck-OrnsteinJITj framework, although it looks for- 
mally quite similar ]q 

The initial conditions for the beam-like partial dis- 
tribution functions are taken to be the same as in the 
nonequilibrium- statistical diffusion approach. Since the 



2 A preliminary version of the bounce back approach had 
been proposed in [12] 
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Fig. 1. Time evolution of pseudorapidity distributions for produced charged particles from minimum-bias d + Au collisions at 
^snn = 200 GeV. Results for five time-steps (p-values, cf. text) are shown. In the left-hand column, the bounce back calculation 
with three sources is displayed. Here the Au- and d-like (grey) distributions reverse their motion in the mean at large times, and 
then tend to separate. Agreement with the data can only be achieved "on the way in" at p=0.30, before bounce back occurs. 
Dash-dotted curves show the midrapidity sources for hadron production. The right-hand column shows the diffusion-model 
calculation where statistical equilibrium would be achieved for large times at the equilibrium value r] eq . The diffusion- model 
evolution reaches agreement with the data much slower, at p=0.54. The minimum x 1 /d.o. /.-values with respect to the PHOBOS 
data [§] are 2.4/48 for diffusion, and 7.4/48 for bounce back. 
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equation is linear, a superposition of the partial distri- 
bution functions using the initial conditions Ri^ilJit = 
0) = S(y ± y m ax) with the maximum rapidity y ma x, and 
R^iVit — 0) = 5(y) yields the exact solution. In the solu- 
tion, the mean values are obtained analytically from the 
moments equations as 

< Vbbi.2{t) >= ±y m ax[l-exp(~t/Ty)]^y max exp(-t/Ty) 

(8) 

for the sources (1) and (2), and 



< Vbb3 (t) >= 



(9) 



for the third source, which rests at in the bounce back 
case during the whole time evolution. It is therefore ob- 
vious that the mean values of the beam-like distribution 
functions start their time evolution at t=0 with values 
Ui2{t = 0) = TUmaxi an d for time to infinity they reach 
yifi{t — > oo) = ±y max . Hence, they change from one ra- 
pidity hemisphere into the other one at a value of the 
bounce back time 



tbb 



= -ln(l/2) = 0.6932 



(10) 



which corresponds to a value of the time evolution pa- 
rameter p=0.5, as can be seen in the middle frame of the 
bounce back case in fig. 1. We use the same t v in sections 
3,4 in order to have a direct comparison of the two sce- 
narios. Of course in the bounce back scenario, there is no 
relaxation towards statistical equilibrium and hence, t v 
should not be considered to be a relaxation time in this 
case. Here it is a time parameter that controls the speed 
of the movement in y- (or 77-)space. 

The variances are as in the diffusion case 



< 2 , 3 W = Dl' 2 ' 6 Ty[l - exp(-2t/r„)]. 



(11) 



The conversion to pseudorapidity is required because 
particle identification is not available: r\ = — ln[tan(0/2)] 
with the scattering angle 9 measured relative to the di- 
rection of the deuteron beam. Hence, particles that move 
in the direction of the gold beam have negative, particles 
that move in the deuteron direction have positive pseu- 
dorapidities. The conversion from y— to r\— space of the 
rapidity density 



dN pdN dN 

~r = = J vn,{m)/{PT))-r- 

drj t, dy dy 



(12) 



is performed through the Jacobian 

(m)/{p T )) = cosh(?7)- 

[l + «rn)/(p T » 2 H 



■ sinh 2 (r])] 



-1/2 



(13) 



We approximate the average mass < m > of produced 
charged hadrons in the central region by a value somewhat 
higher than the pion mass m v , and use a mean transverse 
momentum < pr > = 0.4 GeV/c such that the relevant 
parameter in the Jacobian becomes < m > / < px >= 
0.45/c. 



Due to the conversion from y— to rj— space, the partial 
distribution functions are different from Gaussians. The 
charged-particle distribution in rapidity space is obtained 
in both the diffusion model, and the bounce back case 
as incoherent superposition of nonequilibrium and local 
equilibrium solutions of (|T|) 



dN ch (y,t = Ti n t) 
dy 



= N lh R i(y,n n t) 



+N 2 h R 2 (y, T int ) + NlMy, Tint) (14) 

with the interaction time r int (total integration time of the 
differential equation). In the present work, the integration 
is stopped at the value of Tint /i~ y that produces the mini- 
mum x 2 with respect to the data and hence, the explicit 
value of Tint is not needed as an input. The resulting values 
for Tint/i~ y are given in table YU together with the widths of 
the central distributions, ana the particle numbers in the 
three sources. 

We show the time evolution in the bounce back case 
together with the fit to the PHOBOS data [5] in the left- 
hand column of fig|l] It is evident that the two beam-like 
distribution functions move towards smaller pseudorapidi- 
ties as time increases, reach agreement with the data at 
p—0.30 before motion reversal in the mean is achieved, 
then move on in pseudorapidity space to reverse their mo- 
tion in the mean at p=0.5, and finally proceed until they 
would separate according to the time evolution given by 
eq.Q. At p — 0.999, a very broad distribution in pseu- 
dorapidity has emerged, which deviates considerably from 
the near-thermal distribution that can be seen in the dif- 
fusion case at the same time step. 

The RHIC data of the asymmetric d+Au system show 
very clearly that motion reversal in the mean is not ob- 
served, since the maximum occurs at negative pseudora- 
pidities in the gold direction. Only a certain fraction of the 
produced particles moves opposite to the respective beam 
direction. Although in the bounce back case, the under- 
lying system of differential equations shows a memory in 
the exit channel (P- and T-like distributions emerge sep- 
arately at large times) the effect is not observed, because 
the measurement occurs before that would happen. 

Comparing to the diffusion-model time evolution on 
the right-hand side of fig.l, it is obvious that agreement 
with the data is reached much faster in the bounce back 
case (p=0.30 as compared top=0.54). The reason is found 
in the respective equations for the mean values (|2]), ([8]). 
To first order in 1 1 Tint, the diffusion result is 

< VdXfi >~ TVmax ± y-maxt/Ty (15) 

whereas in the bounce back case 

^ Vbbl,2 >~ -$-ymax 

±2y 

max t/Ty (16) 

and hence, the drift in the bounce back case for fixed t/r y 
is much stronger than in the diffusion case. 

Because the underlying partial differential equation 
conserves the norm, the particle number in each distribu- 
tion function does not change in our schematic model once 
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Fig. 2. Calculated pseudorapidity distributions of produced charged particles in d + Au collisions at ^smn = 200 GeV for 
five different centralities. Central collisions are shown in the top frames, peripheral at the bottom. In the left-hand column, 
results of the bounce back calculation with three sources are shown. The right-hand column gives results of the diffusion-model 
calculation where statistical equilibrium would be achieved for large times at the equilibrium value r/ eq . The diffusion- model 
evolution provides a more detailed agreement with the PHOBOS data [10] . 
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the particles are created. The ^-minimization yields a dif- 
ferent partition of the produced charged hadrons among 
the three sources as compared to the diffusion case, table 
[l] In particular, the percentage of particles in the midra- 
pidity source is much smaller as in the diffusion case, only 

3 % instead of 24 %. As is evident from fig(T|and table [TJ 
the diffusion case yields better agreement with the data. 
However, both the bounce back and the diffusion evolution 
still appear to be acceptable models, and a clear distinc- 
tion between the two becomes apparent only in the sub- 
sequent undetected time evolution. Hence, it seems not 
straightforward so far to conclude from this analysis that 
the system is observed to be on its way to statistical equi- 
librium. To prove that this is likely to be the correct phys- 
ical interpretation, as outlined in [1.2,3,4,5], one indeed 
needs additional physical information. 

A possibility to obtain this information is the investi- 
gation of the centrality dependence as performed in [4 , 5 
for the RDM-case. Such an impact-parameter dependent 
analysis is very sensitive to the details of the model. We 
have performed this analysis for the bounce back case, 
and repeated it for the RDM-case with the same initial 
condition for the central source. We use the same Glauber 
values for the average number of participants in each cen- 
trality bin as in [rUlHIS]. 

The results as presented in fig(2]show that the RDM 
gives consistently better agreement with the data for all 
centralities, in particular for more central collisions where 
the diffusion approach is expected to function best because 
of the larger number of participants. The x 2 — values of 
the bounce back calculation vary from 12.3 for — 20% 
(central) to 4.9 for 80—100% (peripheral), as compared to 

1.5 (central) and 4.6 (peripheral) in the diffusion case. To 
obtain better results in the bounce back case, deviations 
from the Glauber values would be required which are not 
realistic. As an example, for < N\ >= 7 and < N2 >= 1 
in minimum bias we obtain x 2 = 2.1, which is comparable 
to the diffusion result with the Glauber values < Ni >= 

6.6 and < N2 >— 1.7, but the value for the deuteron 
participants is unrealistic. 

The interpretation provided by the Relativistic Dif- 
fusion Model is also favoured by the experimental fact 
that the distributions of produced particles in transverse 
momentum space are very close to thermal distributions. 
Hence, it seems natural to infer that the system tends to 
approach thermal equilibrium, and that this tendency is 
seen in the longitudinal variables as well. 

We have confined this investigation to a comparison 
with the PHOBOS data because these provide a large 
number of data points with high accuracy, giving better 
constraints of the free parameters than the corresponding 
BRAHMS Q3], and STAR results [H]. 

4 Conclusion 

To conclude, we have compared two schematic models for 
charged-hadron production in relativistic d+Au collisions 
at ^/s~nn — 200 GeV for minimum bias, and depending 
on centrality. The nonequilibrium-statistical Relativistic 



Diffusion Model (RDM) describes the gradual approach of 
the system towards statistical equilibrium. The rapidity 
distribution functions become symmetric at large times 
with respect to the equilibrium value of the rapidity, which 
is obtained from energy- and momentum conservation in 
the system of participants. 

Since strong interaction stops before the system reaches 
equilibrium, the pseudorapidity distribution for minimum- 
bias collisions remains rather asymmetric. In this approach, 
all the details of the experimental distribution function 
such as the different slopes in the respective gold and 
deuteron directions are precisely reproduced for all cen- 
tralities with good x 2_v alues. The midrapidity source con- 
tains 24% of the produced charged hadrons. 

In the bounce back approach, the underlying differen- 
tial equation would lead to a re-separation of the partial 
beam-like distribution functions in 77— space, rather than 
to equilibrium. The midrapidity source is very small in 
size (3%). It remains centered at rest, < 773 (i) >= 0. 

When comparing the bounce back time evolution to 
the available PHOBOS-data [5], it turns out that satis- 
factory agreement can only be achieved "on the way in" 
before motion reversal, because the experimental distri- 
bution peaks in the negative pseudorapidity hemisphere, 
fig[T] The data are reached faster during the time evolu- 
tion as compared to the diffusion case because the drift 
is stronger. The system is quite transparent, and it is not 
possible to actually observe motion reversal (bounce back) 
in the mean of the produced-hadron distributions at RHIC 
energies. Most likely, this is true at LHC energies as well, 
since the transparency is expected to be even more pro- 
nounced at higher energies [?. . 

This result is particularly relevant for the interpreta- 
tion of heavy symmetric systems like Au+Au or Pb+Pb, 
which are expected to behave similarly, but where the two 
possibilities - agreement with the data "on the way in" 
(partial transparency), or after motion reversal (bounce 
back) - can not be distinguished easily due to the symme- 
try of the system. From the analogy with our investigation 
of an asymmetric system, however, it appears certain that 
partial transparency is the correct physical description at 
RHIC energies and above. For symmetric systems with 
Deq — 0, the two descriptions outlined in this work are 
difficult to distinguish when compared to data. 

For asymmetric systems, however, the differences of 
the two models are quite pronounced, in particular when 
the centrality dependence is considered. The RDM ap- 
proach is more precise, and we expect that it yields con- 
verging x 2 — minimizations with respect to the experimen- 
tal pseudorapidity distributions for every relativistic heavy- 
ion collision also at other energies, for example, Si+Al 
at AGS or S+Au at SPS energies. The distribution of 
produced particles among the three sources, however, will 
vary substantially depending on incident energy, size and 
asymmetry of the system, and centrality. The third source 
vanishes at sufficiently low energy, for example, at ^/smn 
= 19.6 GeV for Au+Au 0. 

In summary, consistency with the data at RHIC ener- 
gies can only be achieved based on partial transparency 
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without motion reversal in the mean, independently of the 
specific model. Comparing bounce back and nonequilibrium- 
statistical approach in detail, the agreement of the diffu- 
sion approach with the data is definitely better. This result 
supports the view that the relativistic many-body system 
is indeed observed to be on its way towards statistical 
equilibrium. 



5 Acknowledgements 

Two of the authors (MB and TM) would like to thank 
RCNP at Osaka University for discussions at meetings. 
The work is supported by DFG under contract No. STA 
509/1-1. 



References 

1. G. Wolschin, Eur. Phys. J. A 5, 85 (1999); Phys. Lett. B 
569, 67 (2003). 

2. G. Wolschin, Europhys. Lett. 74, 29 (2006). 

3. M. Biyajima, M. He, M. Kaneyama, T. Mizoguchi and N. 
Suzuki, Prog. Theor. Phys. Suppl. 153, 344 (2004); M. Biya- 
jima, M. Ide, T. Mizoguchi and N. Suzuki, Prog. Theor. 
Phys. 108, 559 (2002). 

4. G. Wolschin, M. Biyajima, T. Mizoguchi and N. Suzuki, 
Phys. Lett. B 633, 38 (2006). 

5. G. Wolschin, M. Biyajima, T. Mizoguchi and N. Suzuki, 
Annalen Phys. 15, 369 (2006). 

6. A. Bialas and W. Czyz, Acta Phys. Pol. B 36, 905 (2005). 

7. R. Kuiper and G. Wolschin, Annalen Phys. 16, 67 (2007); 
EPL 78, 22001 (2007). 

8. P.J. Siemens and J.O. Rasmussen, Phys. Rev. Lett. 42, 880 
(1979). 

9. B.B. Back et al., PHOBOS collaboration, Phys. Rev. Lett. 

93, 082301 (2004). 

10. B.B. Back et al., PHOBOS collaboration, Phys. Rev. C 
72, 031901 (2005). 

11. G.E. Uhlenbeck and L.S. Ornstein, Phys. Rev. 36, 823 
(1930). 

12. M. Biyajima et al., preprint |micl-th/0410094y 2 (2004). 

13. I. Arsene et al., BRAHMS collaboration, Phys. Rev. Lett. 

94, 032301 (2005). 

14. B.I. Abelev et al., STAR collaboration, preprint nucl- 
|ex/07030T6| (2007). 



